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Abstract 

Two-pion transitions in charmonium and bottomonium are considered with the 
leading relativistic effects taken into account. The contribution of the chromo- magnetic 
interaction of the charmed quarks to the amplitude of the decay tp{2S) — > ivKJ/ip is 
estimated from the available data. This contribution is enhanced by a factor of three 
in the decay r]c{2S) tttttJc and should produce a noticeable modification of the 
rate and the spectrum in the latter decay. It is argued that the peculiar observed 
spectrum in the decay T(3S') vrvrT arises due to a dynamical suppression of the 
leading nonrelativistic quarkonium amplitude and thus enhanced prominence of the 
relativistic terms. Also discussed are the effects of the final state interaction between 
the pious. 



1 Introduction 



Hadronic transitions between states of the available heavy quarkonia, charmonium and bot- 
tomonium, present a very interesting case study in dynamics of both the heavy quarks and 
the light mesons emitted in the transitions. In particular, the observed properties of such 
processes with the emission of two pions are understood within the chiral low-energy dynam- 
ics of pions, starting with the earliest observation ^l^ of the decay ^'(25') n^n^J/ip and 
the theoretical analyses of the data [21 El- Furthermore, the QCD picture of such decays is 
that the heavy quarkonium transition generates soft gluonic field which then produces the 
light mesons. The heavy quarkonium can be considered as a compact and nonrelativistic 
object in its interaction with a soft gluonic field, which justifies the use of the multipole 
expansion in QCD in analysing these processes jl]. In this approach the amplitude of the 
process factorizes into the heavy quarkonium part, i.e. the transition between the levels as 
a source of the field, and the light meson part, which describes the creation of the mesons 
by the field operators. The former part depends on the dynamics of the quarkonium, while 
the latter one can be understood in some detail by combining [Sj the methods of chiral dy- 
namics and the general low-energy relations in QCD, in particular those for the anomalies 
in the trace of the energy- momentum tensor and in the singlet axial current. Each of these 
two factors in the amplitudes of the discussed hadronic transitions can be used, in certain 
extent, independently. For instance, the decays ip{2S) -n^n^ J/ip and ipi^S) rjj/ip (as 
well as T{2S) tt^tt~T and T{2S) r^T) have essentially the same heavy quarkonium 
part, so that the ratio of the decay rates is fully determined by the low-energy meson am- 
plitudes jni Ej- On the other hand the charmonium amplitude for the double interaction 
with the chromo-electric field (the chromo-polarizability) , which dominates the quarkonium 
factor in the decay 'ip{2S) — > -k^-k^ J/ip can be used [7j for considering the low-energy limit 
of the conversion of ip{2S) into J/ip on a. nucleon: ip{2S) + N J/ip + N, where the soft 
gluonic field is that inside the nucleon. The same charmonium transition amplitude can be 
used as a benchmark for the diagonal chromo-polarizability of J/ip, which determines the 
low-energy elastic scattering of J/ip on a nucleon [HI EI- For these reasons a more detailed 
understanding of both the heavy quarkonium amplitudes and of those determining the gluon 
conversion into light mesons is desirable. 

Presently a quite detailed experimental data are available on the dominant hadronic 
transitions in charmonium 'ip{2S) 71^71^ J/ip, and in bottomonium pUjE]: T(25') —>■ 
TT+TT^T and T(3S') 7i^7i^T{lS,2S). Numerous other, less visible, transitions have been 
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observed, e.g. T{1D) — > tt+tt^T [12j. Thus it is not unreasonable to expect, given the 
capabihties of the modern e+e~ experiments at the charm and bottom thresholds, that a 
significantly improved data will become available. Such data may enable a study of finer 
effects in the amplitudes of the hadronic transitions beyond the leading ones thus far con- 
sidered in the literature. 

The purpose of the present paper is to analyze sub-leading relativistic effects in the am- 
plitudes of the two-pion transitions in charmonium and bottomonium, and possible ways of 
studying such effects from the data. For the transitions between the ^5*1 states of quarko- 
nium, dominated by the second order in the leading El term of the multipole expansion, 
these corrections arise from the second order in the Ml interaction. It will be argued that in 
the charmonium transition ip{2S) — * t^^tx^ J /ip the correction adds coherently to the leading 
term as an 0(10%) term. Being of a rather moderate value in this particular decay, the 
correction is enhanced by a factor of 3 in the transition ric{2S) tt^tt^tIc with a poten- 
tially very noticeable effect in the spectrum of the dipion invariant mass. Furthermore, this 
correction modifies the estimate of the charmonium chromo-polarizability with implications 
for charmonium scattering on nucleons. It will be further argued that, although generally 
the relativistic effects in bottomonium are quite small, in the particular case of the decay 
T{3S) 7r+7r~T it appears that the leading nonrelativistic term it strongly suppressed 
dynamically, and including the relativistic effects of the Ml interaction and of the ^Di Si 
mixing can possibly resolve the long standing puzzle of the peculiar observed spectrum of 
the dipion invariant masses in this decay, which spectrum by far does not conform with the 
predictions from the leading chromo-electric interaction. 

Also some consideration will be given to the significance of the final state interaction (FSI) 
between the pions in the two-pion transitions. The FSI effects have been discussed ever since 
the earliest analyses |2] , and are generally believed to be moderate in the kinematical region 
of the decays ip{2S) T^^T^'J/tp and T(2S') tt+tt^T. However once the amplitudes 
of the transitions are studied in finer detail, one will have to better quantify these effects. 
Although at present the FSI effects cannot be fully analyzed, their behavior at low invariant 
mass of the two pions can be estimated within the chiral expansion. For this purpose we give 
here a systematic derivation of the amplitudes for two pion production by gluonic operators 
in the leading chiral limit, including the terms with the pion mass, and consider the first 
sub-leading term due to iteration of the chiral amplitudes. At the invariant mass of the 
two pions near the threshold these corrections are quite small - only few percent. Their 
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continuation to higher values of the invariant mass can be studied experimentally in a rather 
straightforward way. Such study, in particular, would quantify the impact of the FSI effects 
on the estimates of the quarkonium amplitudes from the observed decay rates. It will be 
also argued that even with the already available data such impact can be estimated as quite 
moderate, as generally expected. ^ 

The material in the present paper is organized as follows. In Section 2 the chromo- 
electric and chromo-magnetic dipole interactions responsible for the two-pion transitions 
are discussed and the parameters determining the relevant quarkonium amplitudes to order 
l(? are introduced. The Section 3 contains a compilative derivation in the leading chiral 
order of the amplitudes for creation of the pion pair by gluonic operators, and in Section 4 
the discussion of the two previous sections is combined in the expressions for the two-pion 
transition amplitudes with the leading relativistic terms included. The latter expressions are 
confronted with the available date on the transitions ip{2S) uTrJ/ip and T{2S) — » vrvrT 
in Section 5 and the significance of the chromo-magnetic interaction in transitions between 
charmonium resonances is evaluated. In Section 6 the FSI effects are considered using 
iteration of the chiral amplitudes and it is also argued on phenomenological grounds that 
these effects are only of a moderate magnitude in the observed two-pion transitions. The 
Section 7 addresses the long-standing puzzle of the unusual spectrum of the dipion invariant 
masses observed in the decay T(3S') vrvrT and it is suggested that this spectrum can 
be explained if the formally leading nonrelativistic quarkonium amplitude is dynamically 
suppressed in this transition, so that the relativistic terms provide a significant contribution. 
Finally, Section 8 contains a summary of the discussion in the present paper. 

2 Multipole expansion for the two-pion transitions and 
relativistic terms 

Considering quarkonium as a compact nonrelativistic object one can apply the multipole 

expansion for the quarkonium interaction with soft gluonic field The leading term 

in this expansion is the chromo-electric dipole interaction with the chromo-electric field E'^ 

^In particular, the present data certainly exclude the recently claimed modification due to FSI of the 
charmonium amplitude by a factor of about 3. 
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described by the Hamiltonian 

HEi = -ler-E''{0) , (1) 

where = — is the difference of the color generators acting on the quark and antiquark 
(e.g. t1 = \°- /2 with A" being the Gell-Mann matrices), and r is the vector for relative 
position of the quark and the antiquark. The convention used throughout this paper is that 
the QCD couphng g is included in normahzation of the gluon field operators, so that e.g. 
the gluon field Lagrangian is written as — (-F^^)^/(4(7^). 

The amplitude of the two-pion transition between heavy quarkonium states 11)2 tt+tt"-?/'! 
(with ■?/'2 and ipi used as a generic notation for the initial and the final quarkonium states) 
can thus be written as ^ 

A{^, ^ vr+TT-^i) = ^-{n+n-\E^E^\0) «Sf , (2) 

(12) 

where alj can be termed, in complete analogy with the atomic properties in electric field, 
as the transitional chromo-polarizability of the quarkonium. In other words, the tp2 —>■ 
transition in the chromo-electric field is described by the effective Hamiltonian 

^eff = "2 "ij^^ E^Ej , (3) 
with the chromo-polarizability given by 

4f = ^(^i|Cr.^r,r|^2) , (4) 

where Q is the Green's function of the heavy quark pair in a color octet state. The latter 
function is not well understood presently, so that an ab initio calculation of the chromo- 
polarizability would be at least highly model dependent. 

Also a simple remark is in order that the discussed decays in quarkonia are governed 
by transitional polarizabilities, i.e. those linking different states of quarkonium. The diago- 
nal chromo-polarizability of quarkonium states, in particular of the charmonium resonances 
ip{2S) and J/ip can also be measured [T3] and is relevant to the problem of scattering of 
these resonances on nuclei [7j. 

Generally aij is a symmetric tensor. Clearly, for transitions between pure S wave quarko- 
nium states this tensor is necessarily proportional to Sij: aij = aoSij. However for the ^5*1 

^For definiteness the processes with emission of a pair of charged pions, tt+tt^ are considered here, since 
the emission of neutral pions tt^tt^ is triviaUy related by the isospin symmetry. 
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states a ^Di — ^ 5*1 mixing generally takes place due to the relativistic effects in the order 
f^/c^ so that a traceless D wave part of the chromo-polarizability is also present and is 
proportional to spin-2 combination of the polarizations, il)2 and ipi of the initial and the final 
states: 

"0 ^ij (V^i ■ ^2) + "2 ^iiV'2j + i'iji'2i - (2/3) 5ij (V^i ■ iJi) ■ (5) 



a 



Numerically, the ^Di Si mixing in charmonium appears to be comparable with the char- 
acteristic value of the relativistic parameter for this system, v'^/c^ ~ 0.2. In particular, 
Rosner jTE] using the e~^e~ decay width of '?/'(3770) and considering this resonance as domi- 
nantly a ^Di state, estimates the angle 9 of the '?/'(3770) — ipi^S) mixing, i.e. the ^Di Si 
mixing for ip[2S), a.s 6 ^ 0.2. It can be noticed however that the mixing does not contribute 
to the spinless part oq of the chromo-polarizability and that in the spin-averaged transition 
rate there is no interference between the S and D wave parts of aij. Thus the effect of the 
^Di — ^ 5*1 mixing in the rate is of order v^/c^ and is thus generally very small. 

It should be also mentioned that the effects of the recoil of the final quarkonium state 
in the decay also arise in the order v'^/c^ and contribute terms quadratic in the momentum 
of the pion pair q to the spin-diagonal part of the tensor polarizability atj. However such 
terms would also be of higher order in the chiral expansion in the pion momenta, which is 
used here. Therefore the recoil terms are suppressed by a product of two small parameters 
and can be neglected within the approximations adopted in this paper. ^ 

Besides the discussed ^Di — ^ 5*1 mixing, the only other effect in the order v'^/c^ in the 
amplitudes of two-pion transitions between quarkonium 5* wave states arises through the 
second order in the Ml interaction with the chromo-magnetic field 5° described by the 
Hamiltonian 

^Mi = -^r(A-S") , (6) 

where A = si — S2 is the difference of the spin operators acting on the quark and the 
antiquark, and M is the heavy quark mass. The contribution of this term to the amplitude 
of two-pion transition between the 5* states of quarkonium can be written as 

Am{^2 TT + TT-^i) = (7) 



^-{n^n-\BtB^\0)aiP 



^ ^ 3 ^ ^ 

5ij (ipl ■1p2) + - ilplii'2j + i'lj1p2i - (2/3) 6ij ipi ■ ^p2) 



■^The recoil effects in tlie quarkonium amplitude were considered in Ref. 17 and estimated to be quite 
small. 
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for the transitions between states, and 

Am{V2 - vr+TT-r^i) = ^(7r+7r-|i?fEf|0)4f (8) 

for the transitions between states. (The notation ri2 and r/i is used here for the states 
of quarkonium in order to emphasize that this relation is specific for the ^5*0 states.) The 

(12) 

amphtude chromo-magnetic analog of the chromo-electric term aij and is formally 

given by the formula 

«!^'^ = 4^(0iir^ri02), (9) 

where (j)2 and (pi are the coordinate parts of the S-wave wave functions of the initial and the 
final states. It is taken into account here that the amplitudes generated by the Ml interaction 
are already suppressed in comparison with those described by Eq.© by the factor v'^/c^. 
Therefore the wave functions of the quarkonium states can be taken in the form where the 
spin and coordinate parts are factorized. 

The spin-2 part of the amplitude in Eq.© contributes to the decay rate only in the 
order w^/c^, similarly to the effect of the ^Di Si mixing. However, the spin-0 part in 
this amplitude as well as that in Eq. (jHJ (the only one present there) does interfere with the 
leading nonrelativistic amplitude in Eq.Q proportional to ao from Eq.(0). Therefore this 
part provides the first relativistic correction, proportional to f^/c^, to the transition rate. 
One can also notice that the effect for the ^5*0 states is three times bigger than for the ^5*1 
states. 



3 Two-pion creation by gluonic operator. Chiral limit. 

As one can see from the previous section, a crucial role in the discussed approach to calcu- 
lating the two-pion transition amplitudes is played by the amplitudes for production of two 
pions by operators quadratic in components of the gluon field strength tensor. Therefore in 
this section we consider the general amplitude of such type: {tt~^ {pi)n~ {p2)\Fl^j^F^^\0) , de- 
scribing the creation of two pions by the local gluonic operator F^^Fx^j. In the leading chiral 
limit the momenta pi and p2 of the pions as well as the pion mass m are to be considered as 
small parameters, and the expression for the amplitude, quadratic in these parameters, can 
be written in the following general form 

- (7r+(pi)7r-(p2)|F;.FAa|0) = [X (pi ■P2) + Y (pI +pI- m^)] {g^xQ.a ~ Q^aQux) + ^ , 

(10) 
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where the structure 



tfiuXa = {PlfiP2\ + PiaP2m) 9ua + iPluP2a + PlaP2u) Q^iX 

- (PlMP2a + PlaP2iJi) 9uX - {PlvP2X + PlXP2u) Qfjia " {Pl " P2) {gfiXQua " gixaQux) (H) 

has zero overall trace: t^^^" = 0, and X, Y, and Z are yet undetermined coefficients. The 
form of the amplitude in Eq. lfTIH) is uniquely determined by the symmetry (with respect to 
the indices) of the operator F^yF\„ and by the Adler zero condition, which requires that the 
amplitude goes to zero when either one of the two pion momenta is set to zero and the other 
one is on the mass shell^. 

The coefficients X and Y in Eq. ()10|) are in fact determined by the anomaly in the 
trace of the energy-momentum tensor O^i, in QCD. Indeed, in the low-energy limit in QCD, 
i.e. in QCD with three light quarks, one finds 

(^' = -^2 C^'^"'^ + E "^.(?"?) > (12) 

q=u,d,s 

where 6 = 9 is the ffist coefficient in the beta function for QCD with three quark flavors. 
The ffist term in Eq. lfT^ represents the conformal anomaly, while the quark mass term 
arises due to the explicit breaking of the scale invariance by quark masses. On the other 
hand, the matrix element of the energy-momentum tensor 9^^ over the pions: 0^u{PiiP2) = 
(7r+(pi)7r~(j92)|6'^jy|0) is fully determined OlIHllini in the quadratic in pi,p2 and m order by 
the conditions of symmetry in fi and u, conservation on the mass shell ((pi+P2)^ dfj.u{pi,P2) = 
at p1 = P2 = m?), normalization {0^u{P:~p) = '^PfiPu at p"^ = m?), and the Adler zero 
condition (6'/,,,(p, 0)|p2=„2 = 0): 

d^uiPi,P2) = [{Pi ■ P2) +pI+pI- m"^] g^lv - Pif,P2u - PiuP2ii ■ (13) 

The equations (fTIHl and (fT!?|l allow for the pion momenta to be off-shell in order to demon- 
strate the Adler zero. However in what follows only the amplitudes with pions on the mass 
shell will be considered, so that it will be henceforth implied that Pi = P2 = ■ In particular 
one finds for the trace of the expression in Eq. ()lH|) 

^;^(pi,P2) =2(pi-p2)+4m2 . (14) 



*The proper index symmetry and the Adler zero condition also automatically ensure that the amplitude 
is C even, i.e. symmetric under the permutation of the pion momenta: pi <->p2. 
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Furthermore, the quark mass term in Eq.()12j). corresponding to the exphcit breaking of the 
chiral symmetry in QCD corresponds to the same symmetry breaking by the pion mass term 
in the pion theory. Thus one finds to the quadratic order in m^: 



(tt+tt I J2 '^g('?'?)|0) = 

q=u,d 



(15) 



while the term with the strange quark makes no contribution to the discussed amphtude. 

Combining the formula in Eq. p2|) for O'^ with the expressions (fT^ and (fT3j) one finds the 
matrix element of the gluonic operator over the pions in the form^ 

327r2 



-{n^p,)n-ip,)\F-F'^>^-^\0) 



2 {Pi ■ P2) + 3 m' 



(16) 



which thus allows to determine the coefficients X and Y in Eq. ()10|) : X = 167r^/(36) and 
Y = 3X/2 = STT^b. 

The coefficient Z of the traceless part in Eq. pO|) cannot be found from the trace relation 
fll2|) . Novikov and Shifman ^H] estimated this coefficient by relating this part to the matrix 
element of the traceless (twist-two) energy-momentum tensor of the gluons only: 6^^^ = 



Zt^A^ = (7r+(pi)7r-(p2)|e|0) 



(17) 



They then assume that the matrix element of the twist-two operator is proportional to the 
traceless part of the phenomenological energy momentum tensor of the pions, 

■1 



(Pl ■ P2) 9f^u - P1mP2z. - PluP2f, 



with the proportionality coefficient interpreted, similarly to the deep inelastic scattering, as 
"the fraction of the pion momentum carried by gluons". They further introduce a related 
parameter k = basPc/iQ'^) conjecture that numerically k ~ 0.15 — 0.20. 

For the purpose of the present consideration the interpretation of the parameter k as 
being related to the pion gluon structure function is not essential and we treat k here as 
a phenomenological and measurable, and actually measured, parameter. Summarizing the 
results so far in this section one can write the expression for the general matrix element (jlO|) 
for on-shell pions as 

.2 r 9 



-(7r+(pi)vr-(p2)|F;,F»JO) 



"3^ 



(q +m) {g^xQv 



(19) 



'^It can be mentioned that this relation, taking into account the pion mass, was used in Ref. (20) . and was 
also derived in a particular chiral model in Refs. |21[I22| . 
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where q = pi + p2 is the total four-momentum of the dipion. 

Few remarks are due regarding effects of higher order in a^. The trace term in Eq. (jl9|) 
receives no renormahzation, provided that the coefficient b is replaced by l3{as)/a1 with 
f5{as) = ba'j. + O(a^) being the full beta function in QCD. This modification however only 
affects the overall normalization of the trace part, and can in fact be absorbed into the 
definition of the heavy quarkonium amplitudes. On the contrary, the relative coefficient of 
the traceless term in Eq. ()19|l . i.e. the parameter k, does depend on the normalization scale, 
which scale is appropriate to be chosen as the characteristic size of the heavy quarkonium jl8j . 
However, given other uncertainties in the analysis of the two-pion transitions, the logarithmic 
variation of k is a small effect. In particular, this effect is likely to be smaller than the 
discussed in this paper relativistic effects in the amplitudes of the two-pion emission. 

The matrix element in Eq. ()19p describes the production of the two pions in two partial 
waves in their center of mass system: the S wave and the D wave. The two waves can 
be measured separately, and also any effects of the final state interaction between pions 
are different in these two orbital states. Therefore it is quite instructive for the subsequent 
discussion to explicitly separate the S and D waves in the matrix element, i.e. to rewrite 
the amplitude ()19p in the form 



Clearly, the trace term in Eq. (jl9|) corresponds to a pure S wave, while the traceless term 
proportional to k, contains both waves. In order to perform explicit partial wave separation 
in t^iuXa it is helpful to introduce the four vector r = pi — P2 describing the relative 
momentum of the two pions, which reduces to a purely spatial vector in the cm. system of 
the pions {{r ■ q) = 0). Then the tensor 



is also purely spatial in the cm. frame and corresponds to pure D wave. Using this tensor 
one can make the following series of replacements for the terms of the generic structure 
PiaP2f3 in the tensor tf,^xa- 




(20) 




(21) 



PlaP2P 



^ Qc^qp - ^ r^rp = - q^qp + — 



4m^ 



) 




(22) 
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Here in the first replacement the cross terms between r and q are dropped since they cancel in 
t^yx^ due to the C symmetry (pi while the gap term in the last transition is dropped, 

since such structure cancels in the traceless tensor t. Using Eq.(j22)) one readily finds from 
the formula (jl9p the expressions for the 5* and D wave amplitudes: 



Si,y\a = ^ {{Q^ + {9fiX9ua - gpLaQuX) (23) 



36 

3 / 2m 



2' 



1 + 



2 V" ' 

and 



Qfiqxgua + quQaQtiX - quqxg,,a - qfiqaQuX {dfiXgua - QfiaQux) 



DfiuXa = i^tiXgua + ^lyaOfiX - ^uXOfia - ^tiaOiyx) ■ (24) 

4 Two-pion transition amplitudes with the relativistic 
corrections 

Using the formulas in the equations 0, and ((Zj) and the expressions p3|l and (j24|l for 
the dipion production amplitudes in the chiral limit, one can readily find the amplitude of 
the transition 1IJ2 — ^ TT^TT'ipi between generic ^Si states of a heavy quarkonium. After a 
straightforward calculation one finds the S wave decay amplitude 



Siij2 ^ TT + TT-^i) = (25) 



"6 



q^ \ 2 



as well as three types of D wave amplitude: one unrelated to the spins of the quarkonium 

states 

Di(^2 VT+TT ^1) = — al ' [1 + xm) — ^^-^^ (^1 ■ i)2) , (26) 

and two amplitudes with the correlation with the polarization of the initial and the final 
resonances 

D2{i^2 - vr+vr-^i) = ^ ^ {x2 + \ Xm) f (l + ^) ^.^.^'^^ (27) 



and 

47r2 / 3 \ 3k 



i^3(^2 - vr+TT-^i) = ^ «f ) (x2 + ^ xm) f ^..^'^^ ■ (2^ 
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In these formulas the following notation is used: P stands for the 4-momentum of the initial 
quarkonium resonance, V'l and ip2 are the polarization 4-vectors for the ^5*1 states, and ip^'^ 
is the spin-2 structure tfj^"" = tp^tp^ + ip^ji^ - (2/3) {ipi ■ {pi'P''/P'^ - g^"'). Finally, xm 
and X2 stand for the ratia 

Xm = , X2 = — (29) 

and encode the relative magnitude of the 0(f^/c^) relativistic effects due to respectively the 
chromo- magnetic interaction (Eq.®) and the ^Di Si mixing. 

The three D waves correspond to different angular correlations. The first one, Di, given 
by Eq.(|26|) corresponds to a D-wave motion in the cm. frame of two pions, which correlates 
with the motion of the cm. system in the laboratory frame, i.e. with the direction of q. 
This D wave arises in the leading nonrelativistic approximation ilSj and is in fact observed 
and measured experimentally jH] for the transition ip{2S) -k^-k^ J /ip. The second D-wave 
amplitude, D2 in Eq.(j2Z|), corresponds to the two pions being in the S wave in their cm. 
system and describes the correlation of the spins of the initial and the final resonances with 
the D-wave motion of the two-pion system as a whole. Finally, the amplitude given 
by Eq. ()28|) corresponds to a D-wave motion of the pions in their cm. frame, which D 
wave is correlated with the spins of the quarkonium states. It can be noted that the two 
latter amplitudes are proportional to a product of two relatively small parameters k and 
0.2 + (3/2) au ~ !(? . Neither nor D3 have yet been observed experimentally, although 
a study [23J of polarization effects in the decay T(2S') tt+tt^T, utilizing a transversal po- 
larization of the DORIS beams qualitatively confirms that these spin- dependent amplitudes 
are quite small. (A discussion can be found in the review [211 •) 

The transitions between ^5*0 states of quarkonium have not been observed yet. One may 
hope however that with a dedicated effort a two-pion transition from the recently found 
rjc{2S) resonance: rjc(2S) n^n^'ric can be observed and studied. Within the approach 
discussed here such transition is closely related to the familiar decay i'i^S) Tr^n^J/ip, 
and in fact can be used for a useful calibration of the total width of ric{2S) |2H1- Clearly, 
on the theoretical side the transitions between ^Sq states are simpler than those between 
the ^5*1 ones since no polarization effects are involved. On the other hand the effect of the 
Ml interaction (Eq.®) is enhanced for the ^Sq states (Eq.©) by a factor of 3, so that the 
relevant transition amplitudes of a generic 772 Tr+Tr^r^i transition are given by 

5(r/2 ^ TT+TT-r^i) = (30) 



11 



47r2 



a, 



(12) 



(1 - 3 Xm) (g' + m^) - (1 + 3 Xm) 1 + 



2m2 



p2 



- Q 



and 



/^i(r?2 



47r2 



a 



(12) 



3k- / PA'P'- 



(31) 



The implications of the enhancement of the relativistic term for the spectrum of the dipion 
invariant masses is discussed in the next section. 



5 Phenomenological analysis 

The most experimentally studied so far transitions of the discussed type are ip{2S) 
71^71^ J/ip and T(25') tt+tt^T. In particular, in the BES experiment [H] both the spectrum 
of the invariant masses of the two-pion system and the angular distribution described by the 
D wave amplitude of the type Di fEq. (pU|) were analyzed, using the formulas equivalent to 
Eq. fl25p and Eq.((2ni) with xm = 0, and the parameter k was determined from the fits. The 
fit to the mass spectrum resulted in the value k = 0.186 ± 0.003 ± 0.006, while the fit to the 
ratio of the Di and S wave amplitude gave k, = 0.210±0.027±0.042. Clearly, the consistency 
of these two values implies that the discussed approach correctly predicts 118^ the ratio of 
the Di wave in terms of the sub-dominant term proportional to k in the S wave amplitude. 
Furthermore, considering the modification of the expressions for these amplitudes due to the 
relativistic parameter xm, it is clear that the fit parameter in the experimental analysis is 
not exactly k but rather 

^'^r=l^^^'^(i+2Ar'*>)^. (32) 

In the estimates in this section we use the value with the smaller error: i^^ff/^^'' = 0.186 ± 
0.003 ± 0.006. 

A similar analysis jTUj of the decay T(2S') tt+tt^T was inconclusive regarding the 
Di wave amplitude, whereas the fit for the dipion invariant mass distribution resulted in 
eff ~ 0-146 ± 0.006. Thus the data display a statistically significant decrease of the 
parameter Kg// in bottomonium in comparison with the similar parameter in charmonium: 
^(^^j/i/.)^^(T^T) _ _ g 27 -|- 0.07. Some decrease of this type of the parameter k was in 
fact predicted ^Hj based on the different characteristic scale for normalization of the gluonic 
operator in Eq.()19|) relevant for the transitions in the two systems. Indeed, the characteristic 
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size of T, rx ~ 1 GeV~ , is somewhat smaller than that of J /tp^ fjl'>p ~ (0.7 GeV)^^. How- 
ever this effect is rather difficult to quantify, given that at both these scales the applicability 
of the perturbation theory in QCD is marginal. Considering also that the difference in the 
size is not large enough, especially in the logarithmic scale, in order to explain the observed 
difference in the values of Kg/ / , one may explore another approach to the observed difference 
by neglecting the variation of the parameter k altogether and ascribing the observed effect 
to different values of the relativistic term xm in these two quarkonium systems. 

Indeed, the relativistic parameter f^/c^ in the lowest states of bottomonium is only about 
one third of that in charmonium, as can be inferred e.g. by comparing the excitation energies 
relative to the mass: (M^(2S) - Mj/^) / {M^^2S) + Mj/^) ^ 3 (Mx(25) - Mt)/(Mt(2S) + Mr). 
Thus one can also expect Xm '^^^"^ ~ 3 Xm ^^"^ assuming that the value of k in Eq.(fT^ is 
the same for the transitions in both quarkonia, one can estimate the relativistic effect of the 
Ml interaction from 

-1^2 xT'^^ - 2x11"'^ ^ 0.27 ± 0.07 . (33) 

Taken at face value, such estimate gives Xm '^'''^^ ~ 0.20±0.05 if Xm is set equal to | Xm '^^^"^ 
and in xit ^'''^^ ~ 0-14 if xiT is neglected altogether. In either case the difference between 
the specific numerical estimates is within the experimental and theoretical uncertainty, while 
the estimated value compares well with the expected magnitude of the relativistic effects in 
charmonium. 

A more definite determination of the parameter Xm '^^^"^ could be enabled by observing 
the distribution in the dipion mass in the decay rjci^S) — >• 'K^Ti~r]c. Indeed, as previously 
mentioned, the effect of the chromo-magnetic Ml interaction is enhanced in this decay by 
a factor of 3. If the value of Xm for the charmonium transitions is in the estimated range 
0.15 — 0.2, the effective parameter i^tf^f"^ ^^e latter decay can amount to 0.3 - 0.45. 
Although the prospect of measuring the Di wave in the decay is likely very remote, the 
effect of such larger value of Kg// should be visible in the more experimentally accessible 
spectrum of the dipion mass, as illustrated in Fig.l. As can be seen from the plots the larger 
value of n^eff"^ results in a significant suppression of the spectrum at low invariant mass, 
due to the zero of the 5* wave amplitude at a higher, than for ip{2S) -k^tt^ J/tp, value of 
g^. Also the total rate contains an overall suppression due to the factor (1 — 3xm)^ vs. the 
factor (1 — XmY for the transition between the ^5*1 states. 



13 



rrin dr_ 

To dm;:;: 




2.5 3 3.5 4 4.5 

mnn/mn 



Figure 1: The spectrum of the invariant masses of the two-pion system described by the 
equations and for the decay ipi^S) rniJ/ip (sohd hne) and for the decay ric{2S) — > 
vrvrr/c with xa4 = 0.15 (dashed) and xm = 0.2 (dash-dot). The rates are normahzed to the 
total rate Fq of the decay i'{2S) nnJ/ip. 

6 Effects of the final state interaction between pions 

So far the amphtudes of the two-pion transitions were considered here in the chiral hmit. 
The formulas in Eq.p9p and in Eqs.()23p and ()24|) are exact in the leading chiral order, 
i.e. as far as the quadratic terms in the pion momenta and mass are concerned. The only 
dynamical modification of these expressions can arise from the previously mentioned QCD 
renormalization effects. In particular these expressions get no corrections due to the final 
state interaction (FSI) between the pions. The latter interaction however can give rise to 
the terms whose expansion starts with the fourth power of momenta and the pion mass, and 
generally can modify the amplitude at momenta of the pions relevant for actual transitions in 
quarkonium. The effects of FSI in chiral treatment of the two-pion transitions in quarkonium 
were a matter of concern ever since the earlier theoretical analyses ^ . The effect in the phases 
of the amplitudes is well known: these phases for the production amplitudes are equal to 
the two-pion scattering phases in the corresponding partial waves: S* = 15*1 exp{i6o), D = 
\D\ exp{i62), where the / = phases for the S wave, 6o, and for the D wave, 62 are quite well 



14 



studied.^ It is also generally estimated both on theoretical and phenomeno logical grounds 
that the FSI corrections are not big (at most 20 - 25%) in the transitions V'(2S') — » vrvr J/-?/; and 
T{2S) —>■ ttttT. (For a discussion see the review [21]-) Some phenomenological arguments in 
favor of such estimate will also be discussed further towards the end of this section, and we 
start with a theoretical estimate of the onset of the higher term in the chiral expansion. 

The interaction of pions at low energy in the D wave is quite weak, so that any modifi- 
cation by FSI of the D wave production amplitude of Eg. ()24|1 can be safely neglected, and 
only the modification of the 5* wave amplitude (j23|l . 5S^u\u, is of interest for present phe- 
nomenology. The imaginary part of the correction at > 4m^ is found from the unitarity 
relation in terms of the isospin 1=0 tttt —>■ tttt scattering amplitude T{q'^) in the 5* wave as 



Im {5S,,^^) = - ^ ^ 5^.A. . (34) 



The amplitude T{q'^) is well known in the chiral limit, i.e. in the quadratic in q and m ap- 
proximation, since the work of Weinberg [22] . In the normalization used here this amplitude 
has the form 

T{,q ) = 7^ , (35) 

J n 

where /^r ~ 130 MeV is the tt"*" fi'^u decay constant. Clearly, the expression in Eg. (j34|l is 
of the fourth power in q and m. 

The real part of SS^^Xa can then be estimated from Eq. ljMj) using the dispersion relation in 
for the amplitude S. In doing so one should set the condition for the subtraction constants 
that this real part does not contain quadratic (and certainly also constant) terms in q and m, 
since these are given by Eq. lf^Hjl . After these subtractions the remaining dispersion integral 
is still logarithmically divergent and contains the well known 'chiral logarithm', depending 
on the ultraviolet cutoff A, which is usually set at A ~ 1 GeV, i.e. the scale where any chiral 
expansion certainly breaks down. Using the equations (|^. and (jH^ one can readily 
estimate the first FSI correction with a 'logarithmic accuracy'. The expression for the full 
5* wave production amplitude then takes the form 

Sf^uXa = + (1 + ^l) i.9^i\9ua - QfiaQux)- (36) 



^It can be mentioned that the analysis [H] of the data on the J/tp ^ n^TT^J/ifi decay does not take 
into account the relative phase between the S and D wave pion production amplitudes. Thus it would be 
interesting to know whether including the phase factor in the angular analysis, produces a significant impact 
on the results. 
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2' 



qf,qx9ua + quqa9^l\ - quq\9f,a - qf^q^gux {9f^\9ua - 9i^a9u\) 



where the correction terms and 62 to Eq.()23p are given as 



^ _ 2(g2)2_7g2^2^3^4 ^^^2 ^_2g2_m2 ^ ^^^^ 



and 



2(g2)2 _9g2^2^g^4 ^2 2g2_m2 / 4m2 
167r2 (9' + 2m2) ^ ^ ' 167r/2 ~ ^ ' 



where the non-logarithmic imaginary part is retained for reference regarding the normal- 
ization. The lower limit under the logarithm is generally a function of both and m^, 
however any difference of this function from the value used in Eqs. flH7|l and ()H8j) is a 
non-logarithmic term, i.e. beyond the accuracy of these equations. Since is the smallest 
of the two parameters in the physical region of pion production, it can be expected that 
using this parameter provides a conservative estimate of the effect of FSI. 

Estimating the corrections in Eq. ()37|) and Eq. ljHHj) . one finds that at the lower end of 
the physical phase space, i.e. near = 4m'^, these terms do not exceed few percent. Thus 
the corrections only weakly modify the normalization of the pion production amplitude near 
the threshold. A theoretical extrapolation to higher values of is problematic, and, most 
likely, one would have to resort to using actual data on the dipion spectra in order to judge 
upon the significance of deviation from the essentially linear in behavior of the amplitude 
described by Eq. (l25|) . A quantitative estimate of the deviation from this behavior has been 
attempted j21j using the data [23] on T(25') tt+tt^T by parametrizing the deviation as a 
factor (1 + q^ /M"^) in the amplitude with M being a parameter. The thus obtained lower 
limit on M is 1 GeV at 90% C.L. 

Another phenomenological argument in favor of a relatively moderate FSI effect in the 
absolute value of the dominant S wave in the decay ip{2S) 7i^7i^J/tp stems from the 
previously mentioned agreement of the observed ^ value of the ratio D/S with the parameter 
K entering the expression for the S wave and extracted from the two-pion invariant mass 
spectrum. Clearly such an agreement would be ruined if there was a significant enhancement 
of the S wave by FSI. 

Furthermore, the observed ratio of the rates of the transitions ip{2S) rjJ/ip and 
ip{2S) 71^71^ J/ ip reasonably agrees with the calculation which neglects any FSI ef- 
fects in the latter decay, while it is clear that there is no such effects in the emission of 
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one-particle, i.e. of rj. The expected accuracy of the theoretical calculation is mostly that of 
applying the nonrelativistic limit to charmonium, i.e. 0(20%) in the amplitude. Thus the 
theoretical (and experimental) uncertainty may allow for a presence of the FSI effects in the 
two-pion system at such level, however large effects of this type, like those recently claimed 
in Ref. [13], are definitely excluded. 

Finally, one can notice in connection with the equations ^I7\ and (j^Hj) for the D2 and D3 
waves involving the polarizations of the quarkonium resonances, that the unknown quarko- 
nium matrix elements all cancel in the ratio D2/-D3. The D2 amplitude is determined by 
the 5* wave of the pions in their cm. system, and the D3 contains the pions in the D wave 
relative to each other. Thus if the ratio of actual amplitudes could be measured from 

angular distributions, e.g. in the decay ip{2S) vr+vr" J/-?/;, this would produce direct data 
on the modification by FSI of the S wave relative to the D wave. 

7 Resolving the T(3S') ttttT puzzle? 

The decay T{3S) ttttT is known to be quite different from the "well behaved" transitions 
ipi'^S) 71^71^ J/ip and T{2S) tt^tt^T in that the spectrum 11, of the dipion invariant 
masses in this decay has two distinct maxima at low and high values of q^. The proposed 
solutions to the puzzle presented by this behavior included presence of a hypothetical exotic 
resonance |2ZIIIZI, breakdown of the multipole expansion [2H1, and unusual FSI effects j29|. 
In this section I exercise an alternative, and somewhat more conventional explanation of the 
observed behavior of the spectrum in terms of the considered here terms of the multipole 
expansion. Namely, one can readily verify that the dipion mass spectrum in the decay 
T{3S) ttttT is reasonably reproduced by the equations (^3)) - (PH)) if the parameters xm 
and X2 for this decay are rather large: of order one, which is probably a result of a dynamical 
suppression of the leading nonrelativistic quarkonium amplitude ao- 

The overall amplitude of the decay T(3S') — * vrvrT is arguably strongly suppressed. In- 
deed, the observed absolute rate of the transition T(35') 7r+7r~T is only about 0.2 of the 
rate of T(2S') tt'^tt'T in spite of a significantly larger available phase space. The estimated 
suppression of the former decay would be even stronger, if its amplitude was "well behaved", 
i.e. given by leading El interaction and thus dominantly proportional to g^. Furthermore, the 
total rate of another similar transition from the same T{3S) resonance, T{3S) — > 7r'7rT(2S'), 
is about 0.6 of the rate of T(3S') — ^ vrvrT, even though the energy released in the former 
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transition is only slightly above the two-pion threshold, and the transition is very strongly 
kinematically suppressed. Thus it appears quite reasonable to conclude that the pure S wave 
part of the quarkonium amplitude in Eq.(0]), ao is small for the discussed transition. This 
can be a result of cancellation in the (appropriately weighed) overlap of the wave functions 
determining this amplitude due to the oscillations of the 35" wave function. Such cancella- 
tion however is not necessarily present in the overlap terms due to the ^Di Si mixing, 
a2 in Eq.(jni), and due to the Ml interaction, aM in Eq-©- The latter terms are naturally 
expected to be of the order 0(f^/c^) in comparison with an unsuppressed leading amplitude, 
e.g. the amplitudes of the "well behaved" transitions: T{2S) —>■ ttttT or T{3S) —>■ 7niT{2S). 
In other words, the ratia xm and X2 (Eq.(j2ni)) can be large due to the small denominator uq. 
A fully quantitative analysis of the decay T{3S) nnT is complicated by that the energy 
released in the transition W = Mx(3S) ~ ~ 895 MeV is rather large for a straightforward 
application of the chiral-limit formulas for the pions. However, for at least a qualitative 
estimate, I neglect here this complication and apply the equations ((23) - to evaluate the 
dipion mass spectrum. The resulting behavior is illustrated in Fig. 2, where the parameters 
are set as xm = 0.7, X2 = 1-0, and k = 0.13. Clearly, the evaluated dipion mass spectrum 
closely resembles the experimentally observed ^1], although no attempt is made here at a 
quantitative fit to experimental data. 

In order to assess whether the values of the relativistic parameters xm and X2 used in 
the plots of Fig. 2 are of the order of the expected relativistic effects in bottomonium, it is 
instructive to compare the corresponding amplitudes a\.j and 0^2 with the amplitude 
of the transition, that appears to be 'normal', namely T(35') —>■ 7r7rT(2S'). In the latter 
transition the released energy is only 332 MeV, and according to the discussed approach it 
should be absolutely dominated by the 5* wave amphtude given by Eq. pH|l . Also any FSI 
effects in this decay should be very small due to the proximity of the pions to the threshold. 
To certain extent this approach can be tested by comparing the rates of the transitions with 
the charged pions and with the neutral, by using respectively the mass of the charged and 
neutral pion in considering each of these decay modes. The effect of the pion mass difference 
is quite essential due to the small available energy, so that after numerical integration of the 
rate calculated from Eq. (j25|) one finds, in place of the isotopic ratio 2, the estimate 

r[T(3^) ^ n-n-T{2S)] ^ 

r[T(35) ^ 7THn{2S)] ■ ' ^ ' 

which is in the agreement with the data B[T{3S) 7r+7r-T(2^)] = (2.8 ± 0.6)% and 
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Figure 2: The spectrum of the invariant masses of the two-pion system as given by the 
equations - for the decay T{3S) vrvrT. Shown in the plot are the total distribution 
(thick line), and the contribution of individual partial waves: S (thin solid line), Di (dashed), 
D2 (dash-dot), and (dotted). 



B[T{3S) 7r°7r°T(25)] = (2.00 ± 0.32)%. ^ The integral over the spectrum of either of the 
modes of transition T(3S') 7r7rT(25') can then be compared with the result of the numerical 
integration over the spectrum of the decay T(3S') — > tt+tt^T produced by the amplitudes in 
the equations ()25|) - ()28|1 which gives the ratio of the rates of the two transitions from T(3S') 
in terms of the ratio of the squares of the corresponding amplitudes ao and can be compared 
with the data. Performing this calculation with the values of the parameters Xm =0.7 
and X2 = 1-0 used in the plots of Fig.2, gives the estimate 

.^^0 ^ ~ n nfi ~ n 04 (do\ 

which certainly falls in the range of naturally expected magnitude of the relativistic terms 

in bottomonium: v^/c^ ~ (^t(2S) — ^^t)/(^^t(25) + ^x) ~ 0.06, and also quantifies the 

dynamical suppression of the leading nonrelativistic amplitude ao in the transition T(35') —>■ 

^In this estimate the value of Kg// used is 0.146. Due to the small energy in the process, the estimate is 
quite insensitive to this particular value as long as te// is small. 
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ttttT. However, given the large energy release in this transition it is rather difficult to quantify 
at present the FSI effect on the estimate (j4(Jj) . In either event, the suggested explanation of 
the observed peculiarity of the decay T(35') —>■ ttttT, heavily relies on the relative prominence 
of the spin-dependent D waves, D2 and D^, in this decay, as can be seen from the plots of 
Fig. 2. The relative contribution of these two waves is the largest near the minimum in the 
invariant mass distribution at 777,7^,7 ~ 0.65 GeV, which is the main experimentally testable 
qualitative prediction of the suggested mechanism. 

It should be mentioned in connection with the suggested presence of the spin- dependent 
D waves that an experimental search for quarkonium polarization effects in T{3S) — >■ ttttT 
has been done ^T]. It was found that the data were consistent at 65% confidence level "with 
the expectation that the daughter T(IS') retains the polarization of the parent T{3S) along 
the beam axis". It thus appears that the test has not been sensitive enough and a more 
detailed study of the bottomonium polarization effects in this transition is needed. 

8 Summary 

In summary. Within the standard approach to two-pion transitions in heavy quarkonium, 
based on the multipole expansion in QCD and chiral dynamics of the pions, the first relativis- 
tic terms of order f^/c^ in the transition amplitudes arise from the effects of the ^Di Si 
mixing in the part determined by the El interaction and from the chromo-magnetic Ml 
interaction. These terms are parametrized by the quantities 02 (Eq-®) and aM (Eq.®)- 
The significance of the chromo-magnetic term in two-pions transitions in charmonium can 
be approximately estimated from the available data on the dipion invariant mass spectrum 
in the decays ipi^S) irnJ/ip and T(2S') vrvrT. The estimated value is 0.15 - 0.2 of the 
dominant nonrelativistic amplitude. The effect of the chromo-magnetic term is enhanced 
in the yet unobserved transition 7]c{2S) — > irnrjc and is expected to significantly distort the 
spectrum in the latter decay and also result in a suppression of its rate, as shown in Fig.l. 
The absolute determination of the leading nonrelativistic amplitude in charmonium, the 
chromo-polarizability, is of interest for other applications, e.g. the charmonium scattering 
on nuclei. Such determination generally suffers from FSI effects of the two pion rescattering. 
These effects are estimated as a next term in the chiral expansion and amount to only few 
percent at a low invariant mass of the two-pion system near the threshold. An extrapolation 
to higher values of can be done using experimental data. With the presently available 
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data there is no indication of a large FSI effect. Furthermore, an agreement of the chiral- 
hmit formulas with the data [0] on the D wave in the decay ip{2S) — ^ irirJ/ilj, as well as 
the agreement with the data of the theoretical prediction [S] for the ratio of the rates of 
the decays ip{2S) rjj/ip and ip{2S) — > T^T^J/ip, suggest that the FSI effects may amount 
to at most a moderate fraction of the amplitude of the transition in charmonium. An ex- 
perimental measurement of the ratio of the D2 wave (Eq.(j2Z))) and the D-^ wave (Eq.(j2HI)) 
in ip{2S) TTTiJ/ip would provide a direct test of the FSI effect. The relativistic terms in 
the two-pion transitions may hold the clue to solving the puzzle of the unusual dipion mass 
spectrum observed in the transition T(3S') — >■ ttttT, if the leading nonrelativistic quarko- 
nium matrix element in this transition is strongly suppressed due to details of bottomonium 
wave functions. Although a detailed quantitative description of the latter decay is not yet 
attainable within the present knowledge, the suggested mechanism necessarily predicts a 
noticeable presence of the polarization-dependent D2 and waves in the amplitude, which 
prediction can be tested experimentally. 
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